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Abstract 

A multi-dimensional extension of the structural default model with 
firms' values driven by diffusion processes with Marshall- Olkin-inspired 
correlation structure is presented. Semi-analytical methods for solving 
the forward calibration problem and backward pricing problem in three 
dimensions are developed. The model is used to analyze bilateral counter- 
party risk for credit default swaps and evaluate the corresponding credit 
and debt value adjustments. 



1 Introduction 

The recent financial crisis has profoundly changed the nature of credit markets 
in general and correlation trading in particular. The focus has shifted from 
complicated products, such as bespoke collateralized debt obligations (CDOs), 
CDOs-Squared, etc., towards simpler products, such as credit indices, collat- 
eralized credit default swaps (CDSs), funded single name credit-linked notes 
(CLNs), CDSs collateralized by a risky bond, etc, for which risks are some- 
what easier to understand and model. However, as recent events have shown, if 
not properly managed, the trading of even these relatively simple products can 



cause big losses. More details can be found in Lipton and Rennie 2011 . 

During the crisis, it has become apparent that proper accounting for coun- 
terparty risk in the valuation of over-the-counter (OTC) derivatives is extremely 
important, especially in view of the fact that some protection sellers, such as 
mono-line insurers and investment banks, have experienced sharply elevated 
default probabilities or even default events, the case of Lehman Brothers be- 
ing the prime example. Counterparty credit risk is the risk that a party to a 
financial contract will default prior to its expiration and will not fulfill all of 
its obligations. In principle, only OTC contracts privately negotiated between 
counterparties are subject to counterparty risk. 



The structural model first introduced by Merton is one of the two standard 
models used for pricing single-name CDSs, the other one being the reduced- 
form model. Extensions of the structural model to the two-dimensional case 
have been proposed by Zhou 2001 , Patras 2006 , Valuzis 2008 , among oth- 



ers, who considered correlated log-normal dynamics for two firms and derived 
analytical formulas for their joint survival probability using the eigenvalue ex- 



pansion technique; see also Lipton 2001 , He et al. 1998 , where an identical 



technique was used in a different context. Two-dimensional structural models 
have been successfully used for the estimation of the credit value adjustment 
(CVA), and the debt value adjustment (DVA) for CDSs (see, e.g 
Sepp||2009|,|Blanchet-Scaillet and Patras| |201l| ). 



Lipton and 



In order to compute the CVA (DVA), one needs to study the joint evolution 
of the assets of the reference name and the protection seller (buyer) , provided 
that the corresponding CDS is viewed from the standpoint of the protection 
buyer. Clearly, in order to calculate the CVA and DVA for a CDS simultaneously 
and consistently, one needs to consider three-dimensional structural models and 
study the joint evolution of the assets of the reference name, the protection 
seller and the protection buyer. This paper extends the results of |Lipton and 



Sepp 2009 by considering correlated log-normal dynamics for three firms and 
computing their joint survival probability. The corresponding problem is solved 
by using the cigenfunction expansion technique combined with the finite element 
method to obtain a semi-analytical expression for the Green's function. Once 
the Green's function is known, both CVA and DVA corrections for a CDS can 
be computed in a consistent manner. The power of the proposed technique is 
illustrated by considering some realistic examples of pricing CDSs sold by risky 
sellers to risky buyers. As might be expected, counterparty credit effects have 
great impact on the value of a CDS contract. 



2 CVA/DVA for CDSs 

A CDS is a contract in which the protection buyer (PB) agrees to pay a periodic 
coupon c to a protection seller (PS) in exchange for a potential cashflow in the 
event of a default of the reference name (RN) of the swap before the maturity of 
the contract T. The value can be naturally decomposed into a coupon leg (CL) 
and a default leg (DL). We denote by t rn the default time of the reference 
name and by Rrn its recovery, and we have (from the protection buyer's point 
of view) : 

CL t = -E [£ T( cD (t,Ti) 1 {Ti < t h» } AT| F t ] , 
DL t = E [ (1 - R RN ) D{t, T RN )l {t<T n. N<T} \ F t ] , 

where are the coupon payment dates and D(t, T) is the price of a zero-coupon 
bond with maturity T. 

In order to simplify the formulas we denote by CF(t,T) the sum of all 
discounted contractual cashflows between t and the maturity T (both coupon 
leg and default leg), and write the value V t of the CDS as: V t = E [CF(t, T)\F t ]. 
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We suppose now that the protection seller can default but consider the pro- 
tection buyer risk free, and denote by Vt the value of the derivative in this 
case: 

V t =E [CF(t,T)l{ T PS >min { T>T RNyy\ F t ] 

+ E [l {T P S<min{TiT H« }} [CF{t,T FS ) + D(t,r ps ) (RpsV+ s + V- PS )} | F t ] , 

where t ps denotes the default time of the protection seller; as usual, V ± — 
±max(0,±V"). The term Credit Value Adjustment (CVA) will refer to the 
additional cost to account for the possibility of the counterparty's default and 
is defined as CVA = V t -V t : 

CVA=(l-R PS )E [l {r P S<min{T!T « N}} C(i,r PS )F+ s | F t ] . (1) 

Similarly we can consider the case where the protection buyer is risky but the 
protection seller is risk free. The term Debt Valuation Adjustment (DVA) rep- 
resents the additional cost to account for one's own default (t pb denotes the 
default time of the protection buyer) : 

DVA = (1 - R PB ) E [t {T PB <min{T>T RH }} D{t, t pb )V; pb | F t ] ■ (2) 

Given recent events, we can no longer suppose that one of the counterparties 
is risk free. The Basel II documentation makes reference to a bilateral counter- 
party risk, in which the default of both counterparties in the derivative contract 
are subject to default risk. One of the advantages of considering the bilateral 
CVA is the symmetry it introduces in pricing: the two counterparties will now 
agree on the price of the derivative (for a detailed discussion on this see for ex- 
ample Brigo and Capponi [2010] ). If r denotes the minimum of the two default 



times: r = min-jV^ , t pb }, then 

V t =E [CF(t,T)l {T>T} + 

+ l {T=T p S<T} (CF(t,r PS ) + D(t,T PS )R PS V+s + D(t,r ps )V- ps ) + 
+ 1 {t=tPB<t} (CF(t,T PB ) + D(t,T PB )R PB V- pB + D(t,T PB )V+ PB )] . 

In the case where both counterparties are considered risky, bilateral CVA is 
the combination of the two adjustments (CVA and DVA): 

CVA= (1 - Rps)E[l {TP s <mm{TiTP B, TRN}} D(t,T PS )V+ PS \ Ft] , (3) 

DVA= (1 - Rp B )E[l {TP B <mm{TtTP s^N }} D(t,T PB )V- pB \ Ft] . (4) 
We emphasize that expressions ([T]), ^ and ([3]), Q are not identical. 

3 Structural model framework 

We assume that the default and counterparty risk can be hedged, so that we 
can work with the risk neutral pricing measure denoted by Q. We also assume 
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a risk- free deterministic rate of return g t . We start with the firm's asset value 



dynamics, which we denote by «/, and assume (similar to the setup in Lipton 



under 



da t 



and Sepp 2009 ) that it is driven by the following jump-diffusion dynamics 

Qt-Ct- A* ft) dt + a t dW t + ((J - 1) dN t , (5) 
a* 

where Q is the dividend rate, W t is a standard Brownian motion, a t is the de- 
terministic volatility, N t is a Poisson process independent of Wt , Xt its intensity, 
j is the jump amplitude and n is the jump compensator. We assume that the 
firm defaults when its asset value becomes less than a fraction of its debt per 
share and that the default barrier of the firm is a deterministic function of time 
given by l t — loE t , where 



E t = exp ( f* (g u - C« - X u n - 



du 



and Iq = RLq. Here R is the recovery on the firm's liabilities and Lq is its total 
debt per share. We consider the firm's equity price per share St and we assume 
that it is given by St = at — h, for t < r and for t > r where r is the default 
time. The solution of the stochastic differential equation ^ can be written as a 
product of a deterministic part and a stochastic exponent a t — loEte (TtXt , where 
the stochastic factor xt has the following dynamics under Q: 

dx t = dW t + — dN t , xo > 0, (6) 

Ot 

with xt representing the "relative distance" of the asset value from the default 
barrier. The default event occurs at the first time r when x T becomes negative, 
so that the default barrier is fixed at zero. 

Introducing jumps in the dynamics of the asset value allows us to calibrate to 
CDS market spreads even for short maturities. In a framework without jumps it 
is well known that a good calibration of the short end of the curve is impossible. 
However, this simpler framework allows for analytical solutions in some cases 
which provide insight into the problem as well as a good benchmark for the more 
general case with jumps. In this paper we therefore focus on the simplified case 
without jumps. 

For the multi-dimensional case we consider that the process for the relative 
distance to default for each of the entities of interest has a similar dynamics to 
equation ([6| but also in the simplified framework without the jump component, 
and we correlate the diffusions in the usual way by assuming d(W£ , W() — Pijdt. 



4 Two dimensional case 

The one dimensional case of the standard CDS where both counterparties are 
non-risky admits well-known analytical solutions. We therefore turn our focus 
directly to the two dimensional problem where we need to model the dynamics of 
the asset values of the reference name and protection seller simultaneously, while 
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considering the protection buyer to be non-risky. We work with the processes 
Xt and yt for the relative distance from the default barrier in time for each 
of the two entities considered. These processes have the following dynamics 
dx t = dWt , dy t = dWjf, where the Brownian motions are correlated with 
correlation p xy . 



4.1 Pricing problem 

The general pricing equation in this framework is given by: 

V t + ^V XX + -Vyy + Pxy V Xy - Q V = 0. (7) 

We consider the following function U(t,x,y) — e e ^ T ~ t ^V(t, x, y) and make a 
change of variable that allows us to eliminate the cross derivatives: 

a(x,y) = x, P(x,y) = -—(p xy x-y), (8) 

Pxy 



where we have used the usual notation p xy = y 1 — p xy . The domain in which 

the equation has to be solved has changed from the positive quadrant to the 
interior of an angle. This angle is characterized by cos(<po) = ~Pxy, so if p xy > 0, 
the angle is blunt. In order to take advantage of the symmetry of the domain, 
we make a second change of variables and pass to polar coordinates, (a, j3) — 
(-r sin (ip -ipo),r cos (ip - tp ))- 



4.2 Green's function 

The Green's function solves the forward problem (where r = T — t) : 



l G +1(3 



0. 



(9) 



0. 



G{Q,r,tp) = ~5(r-r')8{<p-<f/) 
r 

G(r,r,0) = G(r,r,(p o ) = G(t, 0,^=0 G(r,r )V >) 



Two possible methods can be applied in order to obtain the solution for the 
Green's function: the eigenvalue expansion method and the method of images. 
The solution using the first method is well known and was first introduced 
in |Zhou| |2001|, |Lipton| [2001] , |He et al.| |1998|. The resulting formula for the 



Green's function is: 

G(T,r,<p\oy,<p') 



2e~ 



OO 

n=l 



IT 
T 



WKif 



(10) 



A solution through the method of images was announced by Lipton in 2008 
at a SIAM meeting, and briefly discussed in Lipton and Sepp 2009 . Here we 
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present an improved version. We first need to find the solution to equation (|9| 
with the same initial condition but with non-periodic boundary conditions: 

G(T,0,ip) = G(r,r,tp) >0 G{T,r,<p) > 0. 

r->oo — >-oo 

Using a Fourier transform technique, and in order to write the final expres- 
sion in a more compact form we introduce the following function / (p, q) where 

p > 0, — oo < q < oo: 

1 f°° p~P( cosh ( 2 9C)-cos(<j)) 

""O-'-nL—en — *■ 

and its extension h (p, q) defined as follows: 

h (p, q) = 2 [s+f (p, + q) + s -f (p, - q)} , 

where s± = sign (tt ± q). Then we can represent G (t, r, tp\0, r' , </?') in the follow- 
ing form (which can be viewed as a direct generalization of the one dimensional 
case): 

G(T,r,tp\0,r',tp) = —e ^ h l—,<p-tff 

We can now apply the method of images and represent the fundamental 
solution in the form 

oo 

G va (r, r, <p\Q, r', <ff) = £ [G (r, r, ^|0, r', + 2r^ ) - G (r, r, <^|0, r', V + 2rup )] 

71 — — OO 

(11) 



The representation given in equation (11) gives, as expected, exactly the 
same results for the Green's function as those obtained through the eigenvalue 
expansion method. 

4.3 Joint survival probability 

We denote by Q(t,x,y) the joint survival probability of issuers x and y to a 
fixed maturity T. This solves the following pricing equation 

Qt ~t~ ~^Qxx ~t~ ^Qyy PxyQxy 

with final condition Q(T,x,y) = 1 and boundary conditions Q(t,x,0) = and 
Q(t, 0, y) = 0. We use the expression for the Green's function obtained through 
the eigenvalue expansion method and we obtain the following expression for the 
survival probability in the new variables: 

, 4 [°° e-^ 2 ^ 1 /rr'\ . (2k + 1) V 
Q(i, r , 93 ) = - / > ^ Pfc+px — sin rdr. 
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4.4 Application to the CVA computation 



We associate the process x t with the protection seller and the process y t with 
the reference name issuer of a CDS. The protection buyer will be considered 
non-risky in this case. The pricing equation for computing the CVA is given by: 



V t + -V XX + ^Vyy + pxyVxy ~ QV = 0, 



(12) 



with the final condition V(T,x,y) = 0. Boundary conditions are 0, except in 
the case when the protection seller defaults first and there will be a shortfall 
equal to a fraction of the outstanding present value of the standard single name 
swap: 



V 



CVA 



(t,0,y) = (l-R PS )V CDS (t,y) 



where Rps is the recovery of the protection seller. In order to solve this problem 
we apply similar changes of function and variables as in section |4.1| and obtain 
a similar pricing equation 



1 



1. 



1 



U t + - U rr + -U r + -^U, 



2 W 



o, 



with the final condition: U(T,r,(p) — and boundary conditions: 

f7(f, 0, y>) = 0, U(t,oo,<fi) = Q, U(t,r,Q)=0, 
f/ CVA (t, r, p ) = e ^~') (1 - R PS ) V CDS (t, p xy r) + . 

The solution for this problem is given by: 



dG(t' -t,r, (p) 



t Jo 



dG(t' - t, r, ip) 
dip 



U(t',r,0)-drdt' 



ip=0 



(p=ip 



U(t' ', r, (po)-drdt r 
r 



Supplying the boundary conditions for the CVA problem we obtain: 



(13) 



V 



CVA 



(t,r',(p')-- 



l-Rps 
2 



// 

/ 



dG(t'-t,r,ip) 



dip 



V=Vo 



-g(t'-t) v CDS^^ xyr - ) + l drdt , 



(14) 

Thus, by using Green's function, we can perform the CVA computation in a 
natural and straightforward way. 



5 Three dimensional case 

For the three dimensional problem we need to model the dynamics of the asset 
values of the reference name, protection seller and protection buyer simultane- 
ously. As shown previously, we consider the default barrier to be a deterministic 
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function of time for all three assets and we work directly with the processes Xt, 
y t and z t which measure the "relative" distance from the default barrier in time 
for each of the three entities considered. These processes have the following dy- 
namics: dx t = dWf , dy t = dWf , dz t = dW^ , where we correlate the Brownian 
motions with correlations p xy , p xzi p yz . 

5.1 Pricing problem 

The general pricing problem in the octant is: 

Vt + ^V XX + ^Vyy + ^V ZZ + P X yV X y + P XZ V XZ + Py Z Vy Z ~ QV = 0. 

We consider the following function U(t,x,y,z) — e B ^ T ~ t ^V{t,x,y)^ and intro- 
duce a change of variables that allows us to eliminate the cross derivatives: 

n I Pxy% + V) _ ((PxyPyz - Pxz) X + {pxyPxz Pyz) V H~ Pxy%) 

Ot x , p ■ - , 'J 



(15) 



Pxy PxyX 

where we use the notation: \ = Jl — p% y — — p^ z + 2p xy p xz p yz . 

With the change of variables, we have also changed the domain in which we 
need to solve the pricing problem. The domain becomes the volume bounded 

by the planes: a = 0, (0,-^0,7) and (a, ft ?f (- Pxz a + p ' v %'~ Pv ' p) ) ■ 
In order to take advantage of the symmetry of the problem we perform a second 
change of variables to spherical coordinates: the axis a = and j3 — is given 
by 8 = 0; the axis a = and 7 = is given by ip — and 9 — n/2, so that 
(a, j3, 7) = (r sin 9 sin ip, r sin 8 cos ip, r cos 9) . 

The range of values for ip is given by: < tp < ipo where 
<po = arccos (— p xy ). As can be observed in figure [l] the possible range of 



values for 9 will depend on ip, so we have < 9 < @(<p). Formulas (16) and 



(17) give a parametric characterization of this boundary of the domain which 



will prove very useful going forward. 



1 - PxyU 



tp (u) = arccos , (16) 

V y/l - 2p xy u> + uj 2 J 




8 



arccos — - 



Pyz PxzPxy ^ (yPxz Pyzpxy) 



\] Pxy (pIz ~ 2 W {Pxy ~ PxzRyz) + ^ 2 p\z) 

In the domain described above the final form of the pricing equation is: 



Ui 



1 92 
r dr'' 



(rU) 



sm 



2 a 



1 d 

sin 9 89 



(sinOUe) 



= 0, 



(17) 



(18) 



with appropriate boundary conditions depending on the payoff. 
5.2 Green's function 

We now concentrate on solving the forward problem for the Green's function in 
spherical coordinates: 



1 



1^ 

r dr 2 



G(0,r,<p,e) = 



(rG) 



1 



1 



2 \sm 2 9 



G u 



1 d 



(sin 9G e 



0, (19) 



-6 (r - r') S (tp - if') S {9 - 9') , 



r 2 sin 9 

G (r, r, 0, 9) = G(r, r, po, 9) — G (r, r, 0) = 0, 
G(t, r, ya, Q(<p)) = G(r, 0, 9) = 0, G(r, r, if, t 



0. 



We aim to build a solution for the Green's function through the eigenval- 
ues expansion method. The first step is to apply the separation of variables 
technique: 

G{T,r,tp,9)=g(T,r)*{!p,e). (20) 



By substituting (20) in ( |19[ ) we obtain an equation where the left hand side 
depends only on r and r and the right hand side depends only on ip and 9 
and hence both sides are equal to some constant value C, which is necessarily 
negative; we use the notation C = —A 2 . 

For function g(r, r) we have the initial condition g(0, r) = ^5 (r — r') , and 
boundary conditions g(r, 0) = and g(r,r) > 0. Function g solves the 

r— too 

following PDE 

1 fl&_ , , A 2 
r dr 2 



9t 



,.2 ( r 5) - 72 9 



which is similar to the equation utilised for the two dimensional case. The 
solution is given by: 



9(r, 



I 



a/A 2 + 1/4 



In order to obtain the Green's function we also need to solve the angular 
part PDE: 



1 



sm 



. 1 u 



+ 



1 d 
sh^9d9 



(sin 9^e) 



-A 



(21) 
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with zero boundary conditions: #(0,0) = V(tp o ,0) = 0) = *(y>, 8(y>)) = 0. 
It is well-known that the spectrum of this problem is discrete and the set of the 
corresponding eigenvectors is complete. 

The two dimensional spherical surface inside the red line shown in figure [T] 
may be mapped directly onto the (<p, 9) plane. This is done in a similar way to 
the method used by cartographers to map the Earth's surface using Mercator's 
projection. The southern boundary of the domain is mapped into a continuous 
curve parametrised by equations dl6|) and (fl7|). The boundary at = is 



degenerate as it corresponds to the north pole on the sphere. Figure [2] and 
[3] show the domain (denoted by fi) projected onto the {<p, 9) plan for sample 
correlation values. 




Ca " 




Figure 2: p xy 
0.5, p yz = 0.3 



0.8, Pxz 



Figure 3: p xy = 0.8, 
-0.65, p yz = -0.45 



Pxz 



In our case we are interested in writing an eigenvalue expansion for the 
Green's function. The weak formulation for our problem is given by 



1 



sm9V e %dn = A 2 / WsmOdSl, 



(22) 



where VP' is a test function that belongs to the same space as "J, in particular it 
is also on the boundary of the domain. 

The first step necessary is to construct a mesh on the domain of interest. 



We construct a triangular mesh following the ideas presented in Persson 2005 



For the actual mesh generation, the algorithm uses an iterative technique. We 
use adaptive grids that are finer along the boundaries. 

Figure [4] shows an example of mesh obtained for sample values of the correla- 
tions. We show the initial mesh and the final mesh obtained after 100 iterations. 



Once the mesh is constructed we solve the eigenvalue problem in matrix 
form and obtain the eigenvalues and corresponding eigenvectors. Figure [5] shows 
sample eigenvectors for a domain where all three correlations are positive. 
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(a) First iteration mesh 



(b) Mesh after 100 iterations 



Figure 4: Adaptive mesh for the domain obtained for p xy — 80%, p xz = 50%, 
pyz = 30%. The mesh is constructed using 1500 points and is finer as we get 
closer to the boundaries. 



Having calculated the eigenvectors and eigenvalues for our problem we can 
write the eigenfunction expansion for '5 (ip, 9), and then for the Green's function 
we obtain the following final formula: 



G (r, r, <p, 0\ r', ip' , 6') = & —p= ]T ( V ) *" ( ^' 6 )' ( 23 ) 

v n— 1 ^ ' 

5.3 Joint survival probability 

Similarly to the two dimensional case, we denote by Q(t,x, y, z) the joint survival 
probability of issuers x, y and z to a fixed maturity T. This solves the following 
pricing equation 

Qt ~^Qxx ~t~ l^Qyy ~t~ ~^Qzz PxyQxy ~r~ PxzQxz ~T" Pyz Qyz = (24) 

with final condition Q(T, x, y, z) = 1 and zero boundary conditions. We proceed 
with a similar change of variables as described in section 5.1 and using the 
expression for the Green's function given in equation ( 23 ) we obtain: 

OO r 2 +r /2 ^ 



r 5 dr. 



5.4 Application to the CVA computation 

We associate the process Xt with the protection seller, the process yt with the 
reference name and z t with the protection buyer. The pricing equation for 
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(a) Eigenvector 1: = 5.2 



(b) Eigenvector 3: A§ = 16.3 
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computing CVA or DVA in the case where all three names are risky is given by: 
V t + -V xx + ^V yy + -V ZZ + Pxy V xy + Pxz V xz + Pyz V yz -gV = 0, (25) 

with the final condition V(T, x, y, z) = and boundary conditions depending 
on the payoff. 

In the case of the CVA calculation, we get a payout if the protection seller 
defaults, the payout is: 

V CVA (i, 0, y, z) = (1 - R PS ) V(t, y)+, (26) 

where V(t,y) + is the positive value of the single name default swap with non- 
risky counterparts at the time of the default of the protection seller. 
Similarly we have the payout for the DVA calculation: 

V DVA (M,2/,0) = (1 - Rpb) V(t,y)-, (27) 

where Rpb is the recovery of the protection buyer and V(t, y)~ is the negative 
value of the single name default swap with non-risky counterparts at the time 
of the default of the protection buyer. For both CVA and DVA calculations the 
boundary conditions are for all other cases. Following the same procedure as 
in section [5.1 1 we obtain the modified pricing equation 

1 d 2 , s 1 / 1 Id, 

=0, (28) 



U * + \ 



■ or z r z V sin 9 sin 9 o9 



with final condition U(T, r, tp, 9) = and boundary conditions except for: 

U CVA (t, r, 0, 9) =e e ^ (1 - R PS ) V(t, p xy r sin 9)+, (29) 

U DVA (t, r, Vl 9 {ip)) =e^ T -*) (1 - R PB ) V (t, r sin 9 (p xy sin tp + p xy cos tp) ) ~ 

(30) 

for the CVA and DVA respectively. 

The final pricing formula for U is given by: 

U (i, r',tp',9') = - \ JJJ sin (tp) U (f, r, tp, 9 {tp)) Gg (f, r, tp, 9 (tp)) dtpdrdt' 

t 

Tocoo 

IfffU (f, r,<p(u),6(w)) G v (f, r, <p(u),6 (w)) a , , 



VJ J sm6(w) K ' 



t 



2 ,/f .In Jn sinf 



+ 1 rrrmiMG^^^ (3l) 



2 ./ t Jo Jo sin( 
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We note that for one of the integrals above we have used the parametric 



representation of the boundary of our domain given by formulas (16 1 and (17 1 



To obtain the precise formulas for the CVA and DVA calculations we now apply 



the boundary conditions in equations (29 1 and (30) respectively: 



,0(0) 



U CVA (t',r,0,6)G v (f,r,O,0) 



d6drdt', 



(32) 



T oo yo 

1 



U OVA (t,r',^9>)=- - I I /sin0(^);7 DVA (t',r,^0(^))G e (i',r,^0M)d^rdt' 

' (uj)dujdrdt' . 



t o o 

T oooo 



1 f f [U™\1>, r, <p(u),6(u)) GJt', r, <p{u>), qI/ 



V J J sm9 (w) 

t 

(33) 

These original formulas provide a new way of consistently computing the CVA 
and DVA. Similar ideas can be used for many other purposes, which will be 
discussed elsewhere. 



6 Numerical results 

In this section we present the results of the CVA and DVA calculations for a 
risky CDS. We compare the breakeven coupon obtained for a standard CDS 
to those obtained when cither the protection buyer or the protection seller are 
risky (using the 2D formulation and results), as well as when both are risky 
(using the 3D formulation and results). When using the 2D formulation and 
considering that either the protection seller or the protection buyer are risky, the 
two parties will not agree on the breakeven coupon of the CDS. This problem 
goes away when using the full three dimensional framework in which both are 
risky and the problem becomes symmetrical. 

We consider three issuers for our example: AIG as a protection seller, GE 
as a reference name of the CDS, and UNICREDIT a protection buyer. We have 
chosen sufficiently risky entities for the protection seller and the protection 
buyer in order to emphasize the effect of the CVA and DVA adjustments on the 
break even coupon. We calibrate the model inputs to the market data from the 
15th of December 2011 (see table [l]). 

The initial value is a measure of the relative distance to default. The volatil- 
ity a has been calibrated such that the 5Y single name CDS spread is matched 
to the market spread (the 5Y point has been chosen at it is usually the most 
liquidly traded contract). 

For the 2D and 3D cases we also need the correlations between the different 
issuers as inputs to our model. These can be calibrated from the prices of the 
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first to default swap contracts if such contracts including the relevant names 
are available on the market. Alternatively we can proxy these correlations by 
assigning a sector to each issuer and then using the sector to sector historically 
estimated correlations. In this section however we aim to show the impact of 
CVA and DVA on the breakeven spread of a CDS, and hence we will use different 
sets of pairwise correlations for the same group of issuers in order to illustrate 
a variety of cases. 

Figure [6] shows the case where the protection seller is highly correlated to 
the reference name. We observe that the spreads are hyper-exponentially flat 
at 0, which is a known problem for models without jumps. However for longer 
maturities we can match model and market prices. We use the calibrated model 
in order to analyse the effects of either the protection seller, or the protection 
buyer, or both being risky. 

If the protection seller is risky, the probability of it not paying the full 
amount due in the case of the default of the reference name is non-zero, and 
hence the protection buyer pays a lower coupon as it takes on that risk as well. 
If the protection buyer is risky then the breakeven coupon moves in the opposite 
direction and the two counterparties no longer agree on the coupon. Typically, 
the latter shift is much smaller than the former. In the case when both are 
considered risky the problem becomes symmetrical, and a mutually agreeable 
coupon can be computed. 

In this example, in the case of a default of the reference name, the protection 
seller is likely to default as well, and hence the shortfall between the contractual 
payout and what will actually be paid can be significant. The break-even coupon 
will adjust accordingly and will be lower than on a standard fully collateraliscd 
CDS as the expectation of the payout is lower from the protection buyer's point 
of view. 

Figure [7] shows the case where the protection buyer is highly anti-correlated 
to the reference name. This is intuitively the case where the DVA is largest as 
it is in the cases where the reference name does not default that the protection 
buyer is more likely to default on its coupon paying obligation. This leaves the 
protection seller with a potential shortfall. 

7 Conclusion 

A 3D extension of the structural default framework where the joint dynamics 
of the firms' values are driven by correlated Brownian motions is presented. A 



Inputs 


AIG 


GE 


UNICREDIT 


Initial value 


0.0359 


0.3035 


0.1199 


a 


2.44% 


10.45% 


6.3% 


Recovery 


50% 


40% 


40% 



Table 1: Input parameters calibrated to market data (15th December 2011) 
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method to obtain a semi-analytical expression for the Green's function using 
the eigenvalue expansion method is developed. It is shown how to apply this in 
order to compute joint survival probabilities of three different companies and 
how to calculate the credit and debit value adjustments for a standard CDS. In 
the 3D case, a fully analytical expression is not available, since the eigenvalues 
and eigenvectors have to be computed using the finite element method. Given 
a triplet of correlations, however, these can be precomputed, which then allows 
efficient computations across a range of initial points, volatilities or other trade- 
related data (coupons, recoveries etc). Concrete examples demonstrate that the 
CVA and DVA for a typical CDS can be very large. 
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